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t^J- , Abstract The common description of the electrical behavior of a nematic 

ON I liquid crystal as an anisotropic dielectric medium with (weak) ohmic conduc- 

tivity is extended to an electrodiffusion model with two active ionic species. 



Under appropriate, but rather general conditions the additional effects can 
qj", lead to a distinctive change of the threshold behavior of the electrohydrody- 

namic instability, namely to travelling patterns instead of static ones. This 
C*~> ■ may explain the experimentally observed phenomena. 
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1. INTRODUCTION 
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ON ■ In the last years electrohydrodynamic convection (EHC) in nematic liquid crystals (NLC) 

O _ 

has become a standard physical phenomenon for investigating pattern formation and the 

OS ! 

paradigma for anisotropic systems fl|, @, \Q Q|. The main advantage of EHC compared 

O , 

m \ to other pattern-forming systems are the short time scales, the possibility to make cells 



with very large aspect ratios and two easily adjustable control parameters (AC voltage 
V e ff and frequency u! e xt), leading to a wealth of instability scenarios || |J]. For a planar 
(homogeneous) alignment of the director at the confining plates (electrodes) , one observes 
at the instability from the unstructured state (primary instability) normal or oblique rolls, 
i.e. roll axis normal or tilted with respect to the equilibrium director (jj. Both patterns 
can be stationary or travelling |5|, ||, [?], ]8|, [J; furthermore the bifurcation can be continuous 
or (very weakly) hysteretic flO|| . 

The scenario appearing above threshold depends on Lo ex t, the thickness d of the cell, the 
material and the temperature. Especially the conductivity (a material parameter easily 
adjustable by doping), the dielectric anisotropy and, as will be shown below, other electrical 
properties of the material are relevant. 

The commonly used theory to describe EHC are the Ericksen-Leslie equations [11, [T^ l 



combined with the quasi-static Maxwell equations, thereafter referred to as the standard 
model . This theory describes quantitatively (in many cases) the stationary features like 
the threshold voltage, the wavenumber of the rolls and the transition from normal to oblique 
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rolls in particular for sufficiently low external frequencies (conductive regime). The main 
discrepancy between the predictions of the standard model and the experiments, however, 
is the observed Hopf bifurcation. The linear analysis leads for a planar configuration always 
to a continuous bifurcation to stationary rolls [l]]. It is even difficult to produce with this 
model a primary Hopf bifurcation by changing rather drastically some material parameters 
[||. The discrepancies are pronounced for thin cells, pure materials (low conductivity) and 
high external frequencies. 

In this work we replace in the electrical part of the standard model the ohmic con- 
ductivity by two dynamically active species of charge carriers (one positive, one negative). 
The constant ohmic conductivity of the standard model becomes a dynamical variable on 
its own with Ohm's law replaced by migration and diffusion parts of the current. Together 
with the charge density p (coupled to the potential of the induced field via the Poisson 
equation) the new model contains two electrical variables. The new physical processes 
incorporated in the theory are ionic migration, diffusion and dissociation-recombination. 
The corresponding new material parameters are mobilities, diffusion constants and a re- 
combination rate constant. The associated time and length scales are, respectively, the 
migration time T m i g needed for an ion to travel through the cell, the (Debye) diffusion 
length \d and the (linear) recombination time r rec . If r rec is much smaller than the other 
time scales the conductivity can be eliminated and one has p as the only electrical dynamic 
variable. If in addition r m j g is large and is small, drift and diffusion can be neglected 
and we are back to ohmic conductivity and the standard model. Comparing T rec oc l/n 
(no = equilibrium density of charge carriers of each species) and T m i g oc d 2 /V e ff with the 
director relaxation time oc d 2 , and Xd oc n with the roll dimensions, it is evident, 
that our model differs from the standard model just in the parameter ranges where the 
latter fails to describe the experiments. 

In Sec. 2 we present the basic equations and assumptions. In Sec. 3 the model is 
applied to the basic state and linearized around this state. Secion 4 contains an approx- 
imate calculation showing a Hopf bifurcation for sufficiently thin cells, clean materials or 
sufficiently high external frequencies (for materials with negative dielectric anisotropy). 
Sec. 5 gives a discussion and a brief comparison with experimental results. 
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2. BASIC EQUATIONS 

It seems reasonable to assume that the conductivity of NLC's is caused by two species of 
ionic charge carriers arising from dissociation of impurities or of the NLC itself |T^]. We 
assume charges dbe (no loss of generality) and overall neutrality / d 3 rn + = J d s rn~ with 
the number densities n + and n - . The electrical current of each species is caused by 
migration in the electrical fields E, diffusion due to carrier-density gradients, and advection 
with the fluid velocity v, J = J + + J~ with = J_^ ig + «Z«| // + d^dv \)A I • ^ n uniaxial- 
anisotropic systems we have 

Imig = e(^ + n + + ^n~)E =: g E, (1) 
Jdiff = ~e(R + Vn + -R-Vn-), (2) 
ladv = e(n + -n~)v=: pv, (3) 

where p is the charge density and a the tensor of the local conductivity. The mobility 
tensors are pfj = p^dij+p^riirij (n is the director field) with ujjj, //^ essentially independent 
of n + ,n~ and E [13]. We assume that the relative anisotropics for the two ionic species 
are the same and thus equal to the relative conductivity anisotropy, cr a /<7j_, as measured 
in the low-field ohmic range, 

f4j = + —niUj) =: ufa'ij. (4) 

<7j_ 

We adopt the Einstein law 

Df = Vt^, (5) 

where Vr = ksT/e ~ 26mV is the thermal potential. Experimentally this relation is 
satisfied within a factor of about 3 IS, 16, |T^ ]. 

The dissociation-recombination process AB ^ A + +B~ with the undissociated molecules 
or impurities AB and the ions A + , B~ [13] is described in homogeneous systems (no spatial 
variation) by the rate equations h = k^riAB — k r n + n~ (kd, k r =dissociation and recom- 
bination constants) and tiab = —n- In the weak-dissociation limit n + , n~ << uab the 
concentration uab of the undissociated molecules is nearly constant in time, and can be 
replaced by its equilibrium value which we write as Kdi ss n e ^ B = Uq, where Kdi ss = kd/k r is 
the dissociation constant and no is the equilibrium value of n + and n~ . This assumption 
is usually valid when Kdi ss « tiq. For MBBA with ionic dopant TBATPB Ref. [13] gives 
Kdiss = 3.4 x 10 21 m~ 3 while equilibrium densities no = a eq /(e(p + + pT)) are of the order 
of 10 19 m -3 for typical mobility and conductivity values (Tables I and II). 
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TABLE I Scaling 



quantity 



scaling unit 



typical value 



lengths 
time 
voltage 
conductivities 

dependent quantities E and p 
orientational elasticities 
dielectric constants 
viscosities 



d/TT 



.(0) 



(0) 



d 2 



7i 



a eq = pen 
^ and 



t/(0) ,,(0) 2 



d 2 " 

^11 + ^33 



eoe± 
7i = "3 - "2 



-i 



^10...100//m) 



- V 



7T 

0-073 

2.7 V 

10- 7 ...10- 9 (Om)- 



15.3 x 10~ 12 N 
4.2 x Kr n As/Vm 
0.11 kg m s _1 



parameter set MBBA I 



TABLE II material parameters related to conduction 



parameter 



material 



value and source 



mobility p 



mobility ratio 7 



mobility anisotropy fr[ = jj~ 
recombination- 
rate constant k r 



dissociation constant K<n ss 



MBBA 

5CB 
MBBA** 

5CB 
MBBA 

5CB 

dielectric 
liquids 

MBBA** 



0.371 1 [13 



0.06.. .0.35 m, [L7|, I 



1 m 
«i @ 



0.33 Hi; 0.5 |1| 



lOs^n^ 1 

10 15 m 3 g -l 



17] 



14] 



3.4 x IQ^mT 6 |3| 



in units of 10 m /(Vs 



Doping TBATPB. 



In the spatially extended case the rate equations are generalized to 

d t n ± ±-V_-J ± = n ± , (6) 

e 

or with (|) - (D and uab = n e J B 
5 i n ± + V • 

In view of coupling Eqs.@ to the director and momentum-balance equations it is conve- 
nient to write Eqs.(0) as a continuity equation for the charge density p = e(n + — n~) and 



vn 



± 



+ fJ.±g!(±E - V T V)n d 



k r (n,Q 



n + n ). 



(7) 
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a balance equation for the conductivity a := a± = e{p\n + + /u7rt~), 

dtp + Y • [pv + q^Eo - V T ^V (d(7)<r + 2/j, Sl (j)p) = 0, (8) 

<9tc + V ■ cru + pg^E{d(^)a + /isi (7)/)) - pVr^V_ (52(7)0" + d(j)psi(j)p) 

(a + M-p)(cr - 



k r n a eq 



(9) 



where we introduced the effective mobility ,u = /i^ + the mobility ratio 7 = p L /p^_ 
together with d( 7 ) = (1 - 7 )/(l + 7), 81(7) = 7/(1 + 1? and s 2 (7) = (1 + 7 2 )/(l + 7) 2 , 
and the equilibrium conductivity a eq = fieriQ. With the Poisson equation 

P = Y • (e e±e'E), 6' = % + —^^ (10) 

and the potential for S (see Sec. 3) all electrical variables are given in terms of the (scalar) 
potential and a. 

We need five boundary conditions (BCs) for the electrical variables in addition to the 
usual no-slip planar BCs n = (1, 0, 0),v = at the confining plates (z = ±d/2), which act 
as electrodes (NLC layer of infinite extent in the x-y plane) . The applied voltage is given 
by 

d/2 

dzE z = V(t). (11) 

-d/2 

Furthermore there are relations between current, electrical field and density for each species 
at the electrodes which can depend in a complicated way on electrochemical processes and 
may be parametrised e.g. for z = d/2 as Jf = af urface E z - Df urface {nf xt - n±). Two 
limiting cases are i) strongly injecting electrodes, i.e. at least one cr sur f ace is very large, 
which leads to E z = (space-charge limiting conditions) and is (in the isotropic-unipolar 
case) adopted e.g. in Refs |l4|, and ii) no transfer of any charge through the electrodes 
(blocking electrodes), 

J+(z = ±d/2) = J~(z = ±d/2) = 0. (12) 

These last BCs do not involve unknown electrochemical processes and will be assumed 
in the rest of this paper. Blocking electrodes imply that the averaged charge density per 
area, Q := limA-+oo\ Ia^ x dy l^d^^P or ' equivalently, E(d/2) — E{—d/2) is constant 
and we can incorporate permanently adsorbed charges (for future purpose) by setting 
Q = Qad + / dzp = or with the Poisson equation ( |l0| ) 



E z (d/2)-E z (-d/2) = -^. (13) 
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Equations (||)and (g) together with the director and momentum balance, and incompress- 
ibility equations of the standard model fl]], the boundary conditions (|il"l),(|i~2|) and the 
charge adsorption (|T3|), determining the conserved quantity E z (d/2) — E z {— d/2) represent 
the bipolar electro diffusion model. It contains in addition to crj_ and a a three more conduc- 
tion parameters pf\_,"j and k r of the NLC at a given temperature (see Table II). Equations 
(||) and (^) are coupled to the fluid velocity by the advection term and to the director 
by cr' and e' . Conversely, the director and momentum-balance equations of the standard 
model are coupled to the electrical fields via the electrical torque oc CQ€ a (ji • E)E on the 
director and the electrical volume force V T el = pE_ + h.o.t. on the fluid. The director and 
momentum-balance equations do not couple directly to a. 

To see the relative magnitude and time scales of the different processes we scale all 
variables and parameters according to Table I. The result is 

Pi{d t + v- V)p = -V-a'Ecr + DV-a , V_\a- 1 d^)a + 2si(-/)p], (14) 



Pi(d t + v-V)a = -aV_- £E[d{^)a + asx^p] + DV ■ g!V[s 2 (-/)(r + d(j)asi(-/)p] 

Y-e'R = p, (16) 



and scaled director and momentum-balance equations (eqs.(52) in |2C] with a slightly dif- 
ferent scaling). The BCs for an applied voltage V(t) = \/2V e jf cosu> ex tt and blocking 
electrodes are 

r_l) 2 dzE z = v(t), 

[E z a-Dd z (a- 1 d{ 1 )a + 2s l {l)p)] z=±1/2 = 0, (17) 
{Dd z a-aE z {d{ 1 )a + 2as 1 { 1 )p)} z=±1/2 = 0, 

where the scaled voltage is V'(t) = V(t)/V c {0) = V2R cos cot, and we wrote (|12|) as J+ ± 
J- = 0. 

Apart from the two external control parameters R and loq the model contains four 
system parameters P\,a,r and D, see Table III. The first three are (Prandtl-number like) 
time-scale ratios relating the time scales of director relaxation, charge relaxation, carrier- 
generation-recombination kinetics and carrier migration, and D is (up to a factor \/2) 
the diffusion length (Debye screening length) Xd = (pVrT q /2) 1 ^' 2 in units of d/ir. The 
diffusion is not an independent parameter since it is coupled to a via the scaled Einstein 
law D = aV T /v} 0) . 
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TABLE III System parameters of the bipolar model 



parameter 



physical process 



typical value* 



proportionality 



UJQ 



Veff 

w 

, , ,-(°> 

UextT d 



(0) — ^oj 



a 



(v {0) ) 



T, 

r = 



(0), 



= k r n T q 



D = 



n An 

l7T ~ - iy 



control parameter 
control parameter 
charge relaxation 

ion migration 

recombination 
diffusion 



> 2.5 

1...PT 1 
0.057 

1.1 

< 1 

0.1 



* parameter set MBBA I with 
a± = lQ- 8 (nm)- 1 ,n = lQ- 9 m 2 /(Vs), d = lO/xm, and 7 = 1. 



(<J±d 2 



v, 



(0) 



(0) 



nod 2 



a± fie 

(nod 2 )' 1 / 2 oc 



The standard theory is recovered when a — a eq remains negligibly small (and diffusion 
can be neglected). This is the Ccis6 wh.611 i) the recombination time Tree — 

(2/c r n ) 1 is 

much smaller than the charge relaxation time and ii) the ion migration time to travel 
through the cell at the applied voltage is much larger than the charge relaxation time and 
the time scale of the external frequency. This essentially means r >> 1 and a « 1, see 
Table III. 



3. BASIC STATE AND LINEAR ANALYSIS 

In contrast to the standard model the electrodiffusion model actually implies a nontrivial 
basic state ao(z,t), 8po(z,t), E_ = (0, 0, Eq(z, t)), n = (1,0,0) and v = 0. For blocking 
electrodes this is obvious because charge conservation together with finite conductivity 
leads to an accumulation of charges near the electrodes. Other BCs lead to a nontrivial 
basic state (for the electrical variables) as well, except ohmic BCs, J z = cr±E z , implicitely 
assumed in the standard model. 

We decompose the fields of the scaled electrodiffusion model into (i) the basic state of 
the standard model, (ii) the difference between the actual basic state and the standard- 
model basic state (denoted with a S), and (iii) linearized fluctuations, 

n = (1, 0, 0) + (0, n y , n z ) + h.o.t., 

a = a {z,t) + a(r,t) = 1 + Sa (z,t) + a(r,t), 

(18) 

E = E (z, t)e z - V0(r, t) = (ir^V'it) + 5E (z, t))e z - V^r, t), 
p = Sp (z,t) + p(r,t) = Spo(z,t) +i(f> + e a E d x n z + h.o.t., 
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where 5po = d z 5E®, and we used (|l6|), the second equation of ([h]) and the definition 
e = —e' i j(n )didj = — (V 2 + e a 9 2 ) to express p in terms of n z and the induced potential cf). 

Inserting the basic state into (|l4|)-(|l6|) and using 8 z Eq = 5po gives the basic-state 
equations in terms of Eq and the deviations from the trivial basic state, Sao, 5po, 

PidtSpo = | -1 + 2Z?si( 7 )5 2 ] 5p + [-E d z + a^Dd^)^] 5a - 5p 5a , (19) 

Pid t Sa = a (r - l)d(j) - as 1 (j)E d z + Dd(-f)as 1 (j)d 2 5p 

- 2r + ad(j)E d z - Ds 2 (-f)d 2 ] 5a 

- rda 2 + a(r-l)d(-f)5a dp + a 2 (r-l)si(-f)5pl. (20) 

In regions where 5po « 1 and 5gq « 1 the volume charge relaxes with the rate t^ 1 
and diffuses with the effective diffusion constant D^fj^ = 2Ds\(^)d 2 /vr 2 T g = A}? D /T q . The 
deviation from the equilibrium conductivity relaxes with the rate 2r/r q = l/r rec , diffuses 
with the effective constant -De// = Ds2('j)d 2 /ir 2 T q = 2s2("f) X'jy / T q and is advected with the 
effective velocity v a = ad(^)EQd/irr q . Both quantities are coupled by further relaxational, 
advective and diffusive terms. For equal mobilities (7 = 1,^(7) = 0) only the diagonal 
relaxation and diffusion constants and the offdiagonal advective terms survive. The non- 
linear terms oc 5cjq5po and oc da 2 in (^) and fl2Q| ) increase the respective relaxations by 
factors 1 + 5ao and 1 + 5ctq/2. The term oc 5pQ contains the Coulomb repulsion [22] and 
a part of the nonlinear recombination. Equations (^) give the BCs for 5po, o"o = 1 + So® 
and Eq. 



Clearly, the trivial basic state is incompatible with the blocking BCs ( |17D because the 
BCs are not homogeneous in 5po and 5ao. For 'ohmic' BCs J z = a eq E z and — J~ = 
e(pf\_ — pj)E z the BCs ( |17| ) are homogeneous with cjq replaced by 5gq and the trivial basic 
state is a solution. For the blocking or other general BCs boundary layers of accumulated 
charge appear, which influence the impedance of the system. These effects will be discussed 
in future work. Here we assume that their influence on the stability of the basic state is 
small (see below). 

The linearised equations of the electrodiffusion model are obtained by inserting the 
decomposition ( |l8| ) into the scaled equations ( |14| ) and (15). The two electrical equations of 



the resulting set of equations for p and a are (primes on the basic-state quantities denote 
z derivatives and a is defined in analogy to e) 
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Pid t (£(j) + e a E d x n z ) 

= [-a a + a' d z - 2D Sl (~f)ae] <j> - [e 8 z + 5p + a~ l Dd(-/)cr] a (21) 

+ {-(JaVoEo + D oT 1 d{~{)cr a aQ + 2si(^)(a a 5p' - ae a E ) } d x n z - PiSp' v z , 

P x d t a = a {-asi(7) [(E e - 5p' )d z + 5p (a + e + 2re)] 

- d( 7 ) [a a - a' d z + Ds^ae - ra e] } 

- [ad(7)(£ <9 z + *Po(l - r)) + L>s 2 ( 7 )<7 + 2m ] a 

+ {0^(7) [-o- (Ja-E , o + Dsi{i){a a Sp'o ~ ve a Eo) + rcr e a ^ ] (22) 

+ a 2 si(7)-E [-^Poo-a - e a (25p (l + r) + £ <9*)] + Ds 2 (7)cr a crd| 9 s n a - Picr' v z , 

with the (homogeneous) BCs 

[E - a- 1 !)^)^ + [a -2Ds 1 ( 1 )d 2 z ]E z = 0, (23) 
[ J D^-a^( 7 )]^-«K 7 )ao + 2a Sl ( 7 )(<5 / 9o + %)]4 = 0. (24) 

Since for director modes d z ~ n/d and D « 1 the diffusive parts can be neglected leading 
to &/E z = —cjq/Eq and (5po + Eq8 z )E z = 0, which can be fulfilled by the BCs of the 
standard model d z E z = p = for 5po = 0. The linearized director and momentum- 
balance equations of the standard model are changed by some terms from the nontrivial 
basic state. They are given (other scaling) in Ref. Q, Eqs. (3.2) with V(t) replaced 
by V(t) + tt5Eq(z, t) and the volume-force term oc qcj) in (3.2e) {iq corresponds to d x ) 
supplemented by a contribution 5p' Q eoq(j). 

For a — ► 0, D — ► and without the nontrivial parts of the basic state one has a — > and 
Eq. (|2~l"l) reduces to (P\dte + a)(f> = —ir~ 1 (Pie a dt + cr a )V'(t)d x n z , the continuity equation 
(3.2a) in §. 

The BCs for 8po, ao = 1 + 5ao and £b are the same as (|i~7|). Neglecting the small 
field deviations &E7o which are of the order of A^/d, Eq ~ 7r — V'(i) is spatially constant 
and the first two terms on the rhs. of ( JT9| ) and ( p0[ ) are linear in the deviations from the 
standard-model-basic state. 



4. ANALYTICAL RESULTS IN A SIMPLE CASE 



In this section we show that the linear equations (21) - ( |24|) can lead to a Hopf bifur- 
cation. We assume blocking BCs and boundary layers with a thickness much smaller 
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than the cell thickness and in the rest of the volume an essentially trivial basic state, i.e. 
we adopt do = l,nEo(t) = V'(t) = y/2Rcosuit and Spo = 0. Further we neglect diffu- 
sion and restrict ourselves to the normal-roll regime by setting (4>(r,t),a(r,t),n z (r,t)) = 



e zqx (4>(z,t), a(z,t),n z (z,t)) and n y = v y 



0. The fluid velocities v x ,v z in the director 



equation for n z are eliminated adiabatically. This is justified because the viscous relax- 
ation time T v i sc = p m d 2 j '(0:4/2) ~ 10~ 5 s is usually much shorter than the other relevant 
time scales j|, 2C]. With these assumptions one obtains from (H), (||) and the director 
and fluid equations from the standard model 



iqn z 



[Pid t e q + a q }0 + E (t)d z a + [P 1 e a ( J B (i) + E (t)d t ) + a a E (t) 
a [aE (t) Sl (-f)e q d z + d(-f){a q - re q )]0+[Pid t + ad(-f)E (t)d z + 2r]a 
+aE (t) [d(j)(a a - re a ) + asi(7)e a £ , (t)9 z ] iqn z 



cr. 



-ir 2 E (t)<P + 



€ f^El{t) 



iqn z 



0, (25) 

0, (26) 
0, (27) 



normal-roll state) and K q 
is given by pT 



d t - A 

J2 + ( 1 

-<9| + K^q 2 . The zero-field director relaxation operator Aq 



where e g = —d 2 + (1 + e a )q 2 , a q = —d 2 + (1 + cr a )g 2 (the operators e and a applied to the 



Ar 



Ka 



K. 



eff 



-(1-4 



-1 



(28) 



For simplicity the z-dependence of the velocity v z is expressed already here in terms of a 
normalized Galerkin mode f v {z) which gives 



J n<f> 



1-L r 

7I-/2 



q 



tt/2 



dz f*{z) 



«4 + «6 / 2 a2\2 , 2/ 2 



'l+ ai )d 2 z 



(29) 



fv(z) (30) 



with the anisotropic shear viscosities ol^% = 1...6 of the NLC. Some small terms oc have 
been omitted. The BCs (p3|), ( |24| ) without diffusion can be combined to give 

0*^ = 0, d z $ = E (t)W. (31) 



The first BC is consistent with the ansatz for the field inhomogeneity of the symmetric- 
conductive mode of the standard model (Type I A of [Q, |2)). Since then d z <f> / 0, the 
second BC implies the conductivity perturbation to be nonzero at the boundaries, i.e the 
lowest-order mode of the conductivity is antisymmetric in z. 

To get analytical results we use the lowest-order Galerkin ansatz for the symmetric- 
conductive director mode, iqn z = coskz,(f) = ((j)^ e %U)t + e~ tU)t ) cos kz, the anti- 
symmetric mode for the conductivity a = sin kz to which the 4> equation ( |25|) couples 
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n(°) and (f>^ 1 ^ via the Eq8 z term, and the first Chandrasekhar function Q for f v (z). The a 
equation (p6|) couples at lowest order also to symmetric conductivity modes oc e zizlwt cos kz 
but they do not satisfy the second BC (|3l|). Inserting this ansatz into (|25|)-(|27|) one obtains 
four coupled equations for n^°\ 4>^~ l \ 4* an d a ^ ■ Assuming for the charge relaxation 
rate r^ 1 >> T~ e ^,T^ one can eliminate adiabatically the induced potentials c^ 1 ) and 
obtains a 2 x 2 system for the symmetric director /antisymmetric conductivity modes, 

a 

a (o) +{d t -\ n (R,uj)}nW =0. 

The numbers e q ,a q ,K q and Ao are the corresponding operators with d 2 — ► —1. The 
normalized external frequency is to' = oj ex tT q e q / 'a q , where e q /a q varies typically from 0.6 
to 0.9. The relaxation rates X a and A n of the uncoupled a mode and of the conductive 
director mode of the standard model are given by 

K = -(, + "%t-), (33) 



[a,-A„(fl,w)]<T«» + 



Ra 2 (a a ^—e a ) 

* <-> q ' 

71 



{1+U) l2 )cr q <j a 



n(°) =0, 

(32) 



q , 



where R c q is the threshold of the standard model. Besides the two control parameters 
R and u/, Eq. ( |3~2"D contains two parameters describing the migration and recombination 
effects, 

2 _ a 2 si(7) _ ^>l7i7r 2 
Pin 2 a a d z 



(0) 



f = £ = (36) 
1 'rec 

The coupling of the director to the conductive mode is proportional to the square of an 
effective mobility parameter a which is a with /i and T q replaced by the geometrical means 
of fj, ± and T q ,Td, respectively. Experimentally a can be influenced by varying cr±d 2 , see 
Table III. 

The growth rate of modes oc e xt in (|32[) is given by 



C 2 = g V ^"Wl e 9 nnH ) ^ 38 ^ 

where C is (apart from a weak g dependence via q c (u')) a fixed constant for a given NLC. 
For MBBA C 2 (q = q c ) = 1.19 with g c = 1.51 for J « 1 and C 2 (q = q c ) = 2.71 with 
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Figure 1: Hopf frequency ujh = 2tt/h, normalized to the director relaxation time, as 
function of the normalized external frequency to' = oj ex tT q € q l a q for the analytical model of 
Sec. 4 with MBBA I parameters and f = T^/T rec = 1. a is the mobility parameter. The 
cutoff frequency corresponds to uj' = 1.4. 

Figure 2: Hopf frequency as function of the recombination rate for MBBA I parameters 
with e a set to zero (to simulate 152) and different values of a. 

q c — > oo for uj' — > uj' cuto ^. Setting e a = to simulate the material 152, see ||, one has 
C 2 = 0.91, q c = 1.21, independent of frequency (in this approximation). 

A Hopf bifurcation is characterized by a nonzero imaginary part of A at threshold, 
Ke\(R c , q c ) = 0, ImA(i? c , q c ) = u>h / with 



, Rr&C 



1 + uj' 2 ^ C z \a R c a qc J 

The Hopf condition, namely, that the root be positive, leads essentially to conditions for 
f/a and a. For typical material and system parameters (see Tables II and III), one usually 
has ae q /{a q C) « 1 so we are left with an upper limit for the recombination rate, 

aR r C . , . 

1 + Ll>' Z 

Since f/a oc (a±d 2 ) 3 / 2 (Table III) this condition gets weaker for thinner and cleaner cells; 
for materials with e a < it is always fulfilled near the cutoff frequency u' cu t ff, where 
R c — > oo. In the special case e a — > 0, relevant for experiments with the nematic material 
152 |], Eq. (H) simplifies to 



- 2 



27rf H (e a = 0) = -A ay (1 - L nn ) - ^— + (1 - L nn )aj (41) 

If the Hopf conditions are well satisfied, the Hopf frequency is essentially 2nfH = — Aoa(l — 
L nn ) 1/2 proportional to aj_ 1/2 d 3 . For MBBA parameters with e a changed to zero (152) this 
is in physical units fH T d ~ 0.4a. 



5. DISCUSSION 

We presented a new theory for electroconvection by replacing the ohmic conductivity of 
the standard model with drift and diffusion of two species of charge carriers. The con- 
ductivity itself becomes a new dynamically active variable, if the recombination time is 
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not small compared to the charge or the director relaxation time. We demonstrated by a 
simple analytical calculation that this new degree of freedom leads generically to a Hopf 
bifurcation, if the ratio a/r of the mobility parameter a and the recombination parameter 
r is large enough. 

Since the ratio a/r is proportional to (o"j_d 2 )~ 3 / 2 , the above Hopf condition is fulfilled for 
thin and pure (small conductivity a±) samples, in qualitative agreement with experiments 
on MBBA, where Hopf is only found for thin || || and clean cells p3|1 , and with experiments 
on 152 || where the conductivity could be changed and Hopf was found only for low 
conductivities. For materials with e a < (e.g. MBBA) the Hopf condition gets weaker 
(and, if it is fullfilled, the Hopf frequency increases) with increasing external frequencies 
Lu ex t- For e a = the Hopf condition is nearly independent of uj ex t- In fact in 152, Hopf is 
observed in sufficiently clean samples for all frequencies ||. 

In Fig.l our results for the (normalized) Hopf frequency as a function of the (normal- 
ized) external frequency are plotted for MBBA for various values of a and r = 1. The 
behavior found in the experiments of ref.|| is similar to our results for a ~ 0.1 (at r = 1), 
whereas the behavior found in ref.[|| is quite similar to our results for a ~ 0.17. Figure 2 
shows the Hopf frequency as a function of r for MBBA with e a = for three values of a. 

Figure 3: Sketch of the spatial distribution of the physical quantities inside the nematic 
layer (conductive regime). The straight arrows indicate the contributions to the current 
carried by each species. The shading illustrats the conductivity mode (dark = a high, light 
= a low) 
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A qualitative understanding of the mechanism leading to the Hopf bifurcation can be 
obtained by noting that the conductivity mode describing rearrangement of the total charge 
carrier density (each species weighed by its mobility) on the time scale of T m ig reduces the 
charge density produced by the standard Carr-Helfrich charge focussing mechanism. This 
antagonistic cross coupling of o and p can already be identified in Eqs.®,(|) or (|J), dH) 
(one may choose e.g. equal mobilities, 7 = 1, so that for D = 0, r = one is left only with 
cross-coupling terms), and carries through to Eqs,(|32|), where it is easily seen to lead to 
an increase of the threshold (note that in Eqs.(|32|) all effects from the standard theory are 
included in the relaxation rate \ n ). Figure 3 gives a sketch of the resulting distribution of 
the physical quantities inside the layer at two instances separated by half of an external 
period. If the time scale of the cross coupling effect is sufficiently long, the system acts 
like a spring and destabilization occurs in an oscillatory manner. The linear modes of the 
system then describe left and right traveling rolls (or "waves"), although in principle one 
could also have the more complicated superposition to standing waves. 

Since we made crude approximations (adiabatic elimination of the charge, trivial basic 
state, lowest-order Galerkin expansion, normal-roll state also for 152) we can expect the 
above predictions to agree only qualitatively with the experiments. Calculations with the 
full basic state and linear equations are planned in the future. Then it should be possible, 
to determine the recombination rate constant and the mobility ratio of the two species 
by a fit of the predicted Hopf frequency to the measured one. To get information about 
the basic state one can compare the predicted cell impedance with measurements. A 
weakly nonlinear analysis of the new model should determine the parameter ranges where 
the predicted Hopf bifurcation is continuous or hysteretic. Finally one should be able to 
compute all parameters of the universal Ginzburg-Landau description valid near threshold. 
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Figure 1: Hopf frequency ujj = iirfn, normalized to the director relaxation time, as 
function of the normalized external frequency u' = uj ex tT q e q / 'a q for the analytical model of 
Sec. 4 with MBBA I parameters and r = T,i/T rec =1. a is the mobility parameter. The 
cutoff frequency corresponds to u' = 1.4. 




Figure 2: Hopf frequency as function of the recombination rate for MBBA I parameters 
with e a set to zero (to simulate 752) and different values of a. 
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Figure 3: Sketch of the spatial distribution of the physical quantities inside the nematic 
layer (conductive regime). The straight arrows indicate the contributions to the current 
carried by each species. The shading illustrats the conductivity mode (dark = a high, light 
= a low) 
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